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i Abstract. We show that the difference between the Ar 
] and Si relative abundance ratio derived from FUSE ab- 
■ sorption spectra and from the H ii regions of I Zw 18 is a 
consequence of the microturbulent analysis applied to the 
absorption spectra. FUSE observations were performed 
with a large entrance aperture which fully covered the 
galaxy. This means that the observed profiles are averaged 
over the full body of I Zw 18, implying that large-scale ve- 
locity fields influence the absorption-line profiles. Taking 
this into account, we show that the absorption spectra are 
consistent with the same metal abundances as those de- 
rived from the Hii regions. It follows that no significant 
ionization correction as suggested by Izotov and collabo- 
rators to describe metal contents in damped Lya systems 
(DLA) is required to model abundances in the neutral gas 
of I Zw 18 (a local DLA system). Using a mesoturbulent 
approach and applying the generalized radiative transfer 
equation to the Ar i A1048 and Si ii A1020 hues observed by 
Vidal-Madjar et al., we found that the profiles may be re- 
produced with log (Ar/Si) ~ -0.8 and iV(Siii) ~ 4 x lO^^ 

_2 

cm . 
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1. Introduction 

The blue compact galaxy (BCG) I Zw 18 (Mrk 116) 
has been an intensively studied object for the last three 
decades since the first spectroscopic observations by 
Zwicky (1966). I Zw 18 shows an intense and recent burst 
of star formation which makes this galaxy an attractive 
target for studies of star formation history. Recent obser- 
vations revealed, for example, two stellar populations in 
I Zw 18: 10-20 Myr red supergiants and 0.1-5 Gyr asymp- 
totic giant branch stars (Ostlin 2000). 

Amongst BCGs, I Zw 18 shows the lowest oxygen 
abundance : the northwest (NW) and the southeast (SE) 
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bright Hii regions yield log (0/H) = -4.83 ± 0.03 and 
-4.82 ± 0.03, respectively (Izotov et al. 1999). Relative 
to solar values, one finds [0/H][] = -1.70 ± 0.08, i.e. 
Z/Zq ~ 1/50. For argon abundances, Izotov et al. ob- 
tained log (Ar/H) = -6.90 ± 0.05 (NW) and a shghtly 
lower ratio for the SE region, log (Ar/H) = -7.23 ± 0.05. 
Sihcon abundances measured by Izotov & Thuan (1999) 
give log (Si/H) = -6.29 ± 0.22 (NW) and log (Si/H) 
= -6.30 ± 0.22 (SE). Thus, these values imply for the 
NW emission patch [Ar/H] = -1.42 ± 0.06 and [Si/H] = 
-1.85 ±0.22. 

The neutral gas properties in I Zw 18 have been probed 
in both the radio (see van Zee et al. 1998 and references 
cited therein) and the UV range (Kunth et al. 1994; Pettini 
& Lipman 1995; Vidal-Madjar et al. 2000). High velocity 
and high spatial resolution radio observations have shown 
that the overall kinematics of the H i gas associated with 
I Zw 18 is very complex and the neutral gas velocity dis- 
persion a equals 13 — 14 km s^^, or b = \/2a = 18 — 20 
km s^^ (van Zee et al. 1998). The Hi column density in 
front of I Zw 18 deduced from the Lya absorption pro- 
file by Kunth et al., A(Hi) = (3.5 ± 0.5) x lO^i cm'^, is 
comparable to the peak Hi surface density found by van 
Zee et al., A^(Hi) ~ 3.0 x 10^^ cm'^. This means that the 
neutral gas in I Zw 18 can be considered as a local damped 
Lya system (DLA) which is similar to high redshift DLAs 
observed in the light of background quasars. 

First measurement of the 0/H abundance in the Hi 
region have indicated a possible discrepancy between the 
metal content in the neutral gas and in the Hil regions 
(Kunth et al. 1994). This, however, was not confirmed 
in later studies by Pettini & Lipman (1995) and by van 
Zee et al. (1998) who have shown that both the neutral 
and ionized gas in I Zw 18 may have the same oxygen 
abundance. 

New observations carried out with the Far Ultravio- 
let Spectroscopic Explorer (FUSE) by Vidal-Madjar et al. 



^ Using the customary definition [X/H] = log[iV(X) /iV(H)] - 
log[A'^(X)/Ai'(H)]0, and solar abundances from Grevesse et al. 
(1996) except for Ar for which the weighted average value 
-5.48 ± 0.04 from Sofia & Jenkins (1998) is adopted. 
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Fig. 1. (a) - Observed normalized intensities (dots with error bars) vs radial velocities from FUSE spectra of An 
(top panel) and Sin (bottom panel) obtained by Vidal-Madjar et al. (2000). Smooth grey lines are the synthetic 
spectra (convolved with the spectrograph function) of the two-component microturbulent model calculated from the 
minimization as described in the text. Smooth dotted lines show the corresponding unconvolved spectra. Components 
at Av ~ 49 km are the Galactic H2 lines [L(4-0)P(1) in top panel and L(7-0)P(4) in bottom panel] from the 
high- velocity cloud at —160 km s^^. (b) - Same as a but for a one-component microturbulent model, (c) - Confidence 
regions in the Xmin ~ ^ plane calculated from the simultaneous fit of the An and Sill lines shown in panels b. The 
parabola vertex corresponds to the best value of b= 17 A km s~^ 



(2000) produce a similar puzzle : the column density ra- 
tio deduced from the AriA1048 and SinA1020 lines in 
the neutral gas, log (Ar/Si) = —1.32, differs significantly 
from that observed in the Hil regions, log (Ar/Si) nw = 
-0.61 ±0.22 and log (Ar/Si)sE = -0.74 ±0.22 (Izotov et 
al. 2000). Izotov et al. suggested that ionization in DLAs 
may affect abundance ratios. To interpret the observa- 
tions, they suggested a model consisting of two regions 
with substantially different metal contents, i.e. implicitly 
recalling an idea of the existence of two media in BCGs 
(a pristine gas, unprocessed since the big bang and a gas 
polluted by nucleosynthetic products - Kunth & Sargent 
1986). In general, the discovery of such a primordial gas 
would be of great importance for cosmology, as noted by 
Kunth et al. (1994). 

In this Letter, we report on the study of the Ar i A1048 
and SiiiA1020 profiles of I Zw 18 pubhshed by Vidal- 
Madjar et al. (2000). We show that these lines can be 
modeled under the assumption that the metal content in 
the neutral gas is the same as in the H 11 regions without 
refering to an ionization correction. 

2. Data analysis and results 

Spectroscopic observations of I Zw 18 with FUSE in the 
range ~ 980 - 1187 A are described in detail by Vidal- 
Madjar et al. (2000). The spectrum was obtained with a 
resolution of about A/AA ^ 10,000 and a signal-to-noisc 
ratio of S/N ~ 10 per resolution element. The large en- 
trance aperture (30" x 30") fully covers the galactic sur- 
face (- 10" X 4"). As noted by Vidal-Madjar et al., FUSE 
produces the average absorption over the full body of the 



galaxy. In this case the analysis of saturated absorption 
lines is not an easy and unambiguous task. The main diffi- 
culty is connected with the line broadening by large-scale 
irregular (stochastic) velocity fields. The infiuence of the 
finite correlation length on the line profile depends on the 
details of the observation. If one considers the line for- 
mation process in the light of a point source, then the 
observed spectrum reflects only one realization of the ve- 
locity field and, hence, large deviations from the expecta- 
tion value of the intensity may occur if the correlation 
length of the velocity field ^ is not very small compared 
to the size of the absorbing region (Levshakov & Kegel 
1997). If, however, the spectrograph aperture covers an 
essential part of the galactic surface, then should rea- 
sonably well correspond to the observations. But also in 
this case the standard Voigt-fitting analysis (based on the 
assumption of microturbulence) may yield misleading re- 
sults (Levshakov & Kegel 1994). Below we demonstrate 
this effect using the published FUSE data. 

2.1. A microturbulent approach 

We begin with the standard Voigt-fitting analysis which 
assumes a completely uncorrelated velocity field, i.e. (. = Q 
{microturbulence). In our example we use the published 
data from Fig. 1 of Vidal-Madjar et al. to illustrate how 
the column density ratio deduced from two lines, Ar i and 
Si II, is affected by the underlying assumptions. From these 
data given in arbitrary units we have calculated normal- 
ized intensities using the plotted synthetic profiles, and 
then added to the normalized data Icr error bars corre- 
sponding to S/N = 10 (dots and error bars in Figs, la. 
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Table 1. Ar and Si abundances of neutral gas in I Zw 18 
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Fig. 2. Curves of growth for AriA1048 and SiiiA1020 

calculated with the fixed values of the velocity dispersion 
a = 13 km s~^ and the kinetic temperature Tkin = 200 K. 
Three types of lines correspond to micro- , meso- , and 
macrotiubulence {I = 0, L/£ = 1.45, L/£ = 0, respec- 
tively). Dotted horizontal lines mark the estimated values 
of the Ar i and Si ii equivalent widths 



lb). Both profiles were centered at v = km s~^ accord- 
ing to the synthetic spectra of Vidal-Madjar et al., and 

the internal uncertainty of our velocity scale calibration 
was estimated to be about ±5 km s~^. 

The oscillator strengths of the AriA1048.2119 A and 
Sin A1020.6989 A lines, /An 0.257 and /si„ = 0.01391, 
were taken from Federman et al. (1992) and from Charro 
& Martin (2000), respectively. 

The absorption lines of Ar i and Si ii are partly blended 
with the Galactic Ha L(4-0)P(1) and L(7-0)P(4) lines, re- 
spectively. To estimate the Ar i and Si ii column densities 
and the 5-parameter, we firstly fitted two component Voigt 
profiles to the observed intensities via minimization. In 
this procedure the theoretical profiles were convolved with 
the instrumental point-spread function which is supposed 
to be a Gaussian with the width of 26 km s^'^. The best fit 
with xLin = 34.25 (M = 35 data points, and z/ = 28 de- 
grees of freedom) is shown in Fig. la by grey solid curves 
which simultaneously mark the data points involved in the 
optimization procedure. Dotted curves in Fig. la show the 
corresponding unconvolved profiles. As seen from this fig- 
ure, the H2 blends do not affect the cores of the Ar i and 
Si II lines much. This is why in the further analysis we used 
only the blue wings and the central parts of these lines. 

Our main results are shown in Figs, lb and Ic. Fig. lb 
presents the best fit with Xmin = 18.36 (M = 19, = 16), 
Nati = 6.83 X 10^3 cni-2, Nsm = 1-11 x lO^^ cm^^, and 
b = 17.4 km s"-"^. Fig. Ic shows the calculated Xmin values 
as a function of b with the la and 2a confidence levels. 
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^ the internal accuracy of these values is ~ 0.1 dex 



Here, for each point we fixed the value of b and minimize 
by varying all other parameters. 

The estimated ratio of log (Ar/Si) = —1.21 is consis- 
tent within the imcertainty interval of 0.11 dex with the 
value —1.32 cited by Izotov et al. (2000). The most likely 
value for b, 17 A ± 3 km s~^, is in an excellent agreement 
with the radio observations of the H i gas in I Zw 18 car- 
ried out by van Zee et al. (1998). The fact that 6hi — b^vi 
implies that the kinetic temperature of the neutral gas in 
I Zw 18 is less than 600 K and that the line broadening is 
caused mainly by turbulent motions. 

To summarize this section we note that the microtur- 
bulent results lead to the puzzle mentioned above which 
encouraged us to try a more general model. 



2.2. A mesoturbulent approach 

21 cm observations with high spatial resolution show that 

the line profile varies with position (van Zee et al. 1998). 
This is a clear indication that large-scale motions deter- 
mine the line profiles and that the microturbulent ap- 
proach is not well founded. We therefore generalize our 
analysis to include the effects of a finite correlation length, 
i.e. 1^0 {mesoturbulence) and to study the limiting case 
L/£ (macroturbulence) . 

The simulation of the mesoturbulent Ar i and Si 11 pro- 
files has been carried out using the simplified model of 
a plane-parallel slab of geometrical size L with homoge- 
neous turbulence and uniform kinetic temperature, Tkin- 
Specifying the An and Sin column densities, the L/£ ra- 
tio, the velocity dispersion a, and the thermal widths for 
each line, we can calculated their average absorption-line 
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profiles employing the generalized radiative transfer equa- 
tion (Levshakov & Kegel 1994, 1997). 

The same ^^-minimization procedure as before was ap- 
plied to the same data set, but now with two fitting pa- 
rameters {iVsiii, log(Ar/Si)}, i.e. v = 17. We fixed a — li 
km s~i and Tki„ = 200 K (i.e. 6 ~ 18 km s^^) and cal- 
culated Xmin ^ given L/l. The results are presented in 
Table 1. It should be emphasized that all Ar i and Si ii pro- 
files corresponding to the listed - very different - solutions 
are identical to those shown in Fig. lb. 

Our numerical results show that accounting for corre- 
lation effects strongly affects the derived column densities 
as well as the abundance ratio Ar/Si. These effects are 
closely related to the changes in the curves of growths 
caused by the finite L/l value. Fig. 2 shows three curves 
of growth for each of the two lines. They correspond to the 
two limiting cases of micro- and macroturbulence (Gail et 
al. 1974), as well as to an intermediate - mesoturbulent 
- case {L/l = 1.45). In general, one may say that the 
column density derived from a measured equivalent width 
increases with increasing correlation length I. Fig. 2 also 
allows a qualitative interpretation of the Ar/Si ratio in 
dependence of L/ i. In the microturbulent limit both lines 
are still close to the linear part of the curve of growth, 
while in the mesoturbulent regime they lie on the flat part. 
Since the slope here is smaller than on the linear part, the 
Ar/Si ratio is larger. In the macroturbulent limit, both 
lines lie on the square root part. In this case the slope 
is again steeper and, more importantly, the two curves of 
growth are well separated from each other (due to the dif- 
ferent damping constants), implying a lower Ar/Si ratio 
than in the mesoturbulent case. Thus, going from micro- 
to macroturbulence we expect the Ar/Si ratio at first to 
rise and, after going through a maximum, to decline again. 
From Fig. 2 we also see that in the macroturbulent limit 
the Ar/Si ratio is lower than in the microturbulent one. 

3. Conclusions 

We have shown that metal abundances derived from the 
absorption profiles of AriA1048 and SiiiA1020, observed 
in the spectrum of the damped Lya system I Zw 18 by 
Vidal-Madjar et al. (2000), become fully consistent with 
those derived from emission-line spectra of the H ii regions 
by Izotov & Thuan (1999) and by Izotov et al. (1999), if 
correlations in the large-scale velocity field are accounted 
for. The adequacy of the obtained mesoturbulent solutions 
is supported by the fact that the derived velocity disper- 
sion of the neutral gas, cr ~ 13 km s""'^, is in an excellent 
agreement with the Hi 21 cm observations, cr = 12 — 14 
km s~^ (van Zee et al. 1998). This implies that the gas in 
I Zw 18 is efficiently mixed. 

To estimate more accurately metal abundances in the 
H I gas within the framework of our model, one needs to 
know the L/ i ratio. This parameter can be found from the 
analysis of saturated and optically thin lines of the same 



ion, since the latter are less affected by the correlation 
effects. For our case, observations of the AriA1066 line 
would be of particular interest since its oscillator strength 
/i066 = l/4/io48- The AriA1066 line was observed in the 
z = 3.4 DLA system toward the quasar QOOOO-2620 and 
its analysis showed a remarkably similar abundance to the 
other a-chain elements O, S, and Si, Z/Zq ~ 1/80 (Molaro 
et al. 2000). The similarity of both DLA systems is also 
supported by the very low molecular hydrogen contents 
found in I Zw 18, /(H2) = 2Nb.JNb.i < lO'^ (Vidal- 
Madjar et al. 2000), and at z = 3.4, /(H2) ~ 4 x 10"^ 
(Levshakov et al. 2000). 

Since Ar can hardly be depleted onto dust grains, but 
can be partially ionized by nearby UV stellar radiation 
with energies hv > 15.76 eV, the relative abundance of 
Ar I is a good indicator of the intensity of the local pho- 
toionizing fiux (Sofia & Jenkins 1998). Our measurements 
rule out the presence of a significant amount of partially 
ionized gas in the damped La system I Zw 18 and, hence, 
metal abundances in the neutral gas do not require ion- 
ization corrections as suggested by Izotov et al. (2000). 

Thus, it is extremely important in future observations 
to investigate low-ion lines to undestand better the kine- 
matic characteristics of the neutral gas bulk motion which 
will enable us to obtain more reliable estimations of mctal- 
licities. 
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